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Abstract 

We study the exact solutions of the cascade three-level atom interacting with a single 
mode classical and quantized field with different initial conditions of the atom. For the 
semiclassical model, it is found that if the atom is initially in the middle level, the time 
dependent populations of the upper and lower levels are always equal. This dynamical 
symmetry exhibited by the classical field is spoiled on quantization of the field mode. To 
reveal this nonclassical effect an Euler matrix formalism is developed to solve the dressed 
states of the cascade Jaynes-Cummings model (JCM). Possible modification of such effect 
on the collapse and revival phenomenon is also discussed by taking the quantized field in 
a coherent state. 
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I. Introduction 



Over the decades, studies of the population inversion of the two, three and multilevel 
system have been proved to be an important tool to understand various fundamental aspects 
of quantum optics [1,2]. Many interesting coherent phenomena are observed if the number of 
involved levels exceeds two. In particular, the three level system exhibits a rich class of coherent 
phenomena such as two photon coherence [3], double resonance process [4], three level super- 
radiance [5], coherent multistep photo- ionization [6], trilevel echoes [7], STIRAP [8], resonance 
fluorescence [9], quantum jump [10], quantum zeno effect [11] etc [12-16]. From these studies, 
it is intuitively clear that the atomic initial conditions of the three level system can generate 
diverse quantum optical effects which are not usually displayed by a two level system [17-20]. 
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The idea of the present investigation is to enunciate the three level system for various initial 
conditions while taking the field mode to be either classical or quantized. In this paper the 
three level system is modelled by the matrices which are spin-one representation of SU(2) 
group. A dresscd-atom approach is developed where the Euler matrix is used to construct the 
dressed states. Wc discuss the time development of the probabilities both for the semiclassical 
model and the cascade JCM for various initial conditions and point out the crucial changes. 
Finally the collapse and revival phenomenon is presented taking the quantized field initially in 
a coherent state. 

The subsequent Sections of the paper arc organized as follows. To put our treatment in 
proper perspective, in Section 11 we have derived the probabilities of three levels taking the field 
as a classical field. The cascade JCM and its solution in the rotating wave approximation (RWA) 
is presented in Section III. In section IV we have numerically analyzed the time dependent 
atomic populations and compared with the semiclassical situation by taking the quantized field 
initially in a number state and in a coherent state. Finally, in conclusion, we highlight the 
outcome of our paper and make some pertinent remarks. 

II. The Semi-classical Model 

The Hamiltonian to describe the semicalssical problem of a cascade three level system 
interacting with a single mode classical field is 



where /'s represent the spin-one representation of SU{2) matrices corresponding to the cascade 
three level system with equal energy gaps(/ia;o) between the states, namely. 
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TvjJi is the interaction energy between the three level system with the classical field mode of 
frequency u in RWA. Let the solution of the Schrodinger equation. 



with Hamiltonian (1) is given by 

iP{t) = C+{t) I + > + Co(t) I > + C_(t) I - >, 
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where C+(t), Co{t) and C_(t) are the time dependent normalized amplitudes with the eigen 
functions given by 
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We now proceed to calculate the probability amplitudes of the three states. Substituting Eq. (4) 
into Eq.(3) and equating the coefficients of | + >, | > and | — > from both sides we obtain 



ihC+{t) = ouoC+{t) + ^exp(-iu;i)Co(t) 
ihCo{t) = ^ex.p{iujt)C+{t) + ^ex.p{-iut)C^{t) 

ihC-{t) = ^exp{iu;t)Co{t) - uJoC_{t), 

where the dot represents the derivative with respect to time. 
Let the solutions of Eqs.(6) are of the following form, 

C+(t) = ^4+ exp(is+i), 
Co{t) = Aoexp{isot), 
C-{t) = A^exp{is-t), 
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where A's are the time-independent constants to be determined. Plucking back Eqs.(7) in 
Eqs.(6) we obtain 



{sq-UJ + LlIo)A+ + -^(^iAq = 

sqAo + -^u;{A+ + A_)^0 

{so + cu - cJo)A- + -^^^iAq = 0. 
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In deriving Eqs.(8), the time independence of the amplitudes A+, A_ and Aq is ensured 
by invoking the conditions s+ — sq — cu and s_ — sq + cu. The solution of Eqs.(8) readily yields 

So = (9a) 

So = -u;of + ujf] (= ±0) (9b) 
and we have three values of s+ and s_, namely, 

s\ — —cu, = fl — LU, s+ = — <^ (10a) 

sl^u, s'i^n + u;, = + (10b) 

Using Eqs.(lO), Eqs.(7) can be written as 

C+{t) ^ A\_ex.p[-iu;t] + Alex.p[i{Q - u;)t] + Alex.p[i{-^1 - u;)t] (11a) 

Co{t) = Al + Al ex.p{int) + Al exp{-int) (lib) 



C-{t) = exp{iujt) + Al exp[i{n + uj)t] + A^_ exp[i(—fl + u!)t]. (11c) 

where As' be the constants to be calculated from the following initial conditions: 

Case-I : Let us consider at t = 0, the atom is in the lower level i.e, C+(0) = 0, 
Co(0) = 0, C_(0) = 1. Using Eqs.(6) and (11), the time dependent probabilities of the three 
levels are given by 

|C+(i) |2 = ^sin^Qi/2, (12a) 

\Co{tW = + (12b) 

|C_(i) p = ^[{ujlsm'^^}t/2 + Q2cosQi)2 + {ou - uJo)^n^sin'^nt]. (12c) 

Case-II : If we choose the atom is initially in the middle level, i.e, C+(0) = 0, Co(0) — 
1, C_(0) = 0, the corresponding probabilities of the levels are given by 

|C+(i)p = ^[4(a;-u;o)^sin^^^t/2 + Q2sin2Qt] 

= \C-{t)\', (13a) 

\Co{t) p = ifc9^2l^sin^QV2 

+^^^^g^\m^flt/2cosnt + cos'^flt. (13b) 

Here we note that, unlike the previous case, the probabilities of the upper and lower levels are 
equal. 

Case-Ill : When the atom is initially in the upper level i.e, C+(0) — 1, Co(0) = 0, 
C_(0) = 0, we obtain the following occupation probabilities in the three levels, 

\C+{t) \'^^:^[{ulsm^nt/2 + n'^cosnty + {u-uJofn'^sm'^nt], (14a) 
\Co{t) \^ ^ ^[4{u-uoysm'^flt/2 + fl'^sm^flt], (14b) 

|C_(i) |2 = ^sin^Ot/2. (14c) 

We note that the probability of the middle level for case-Ill is precisely identical to that 
of case-I while those of the upper and lower levels are interchanged. 

III. Cascade Jaynes-Cummings Model 

Here we consider the cascade three level system interacting with a single mode quantized 
field. The cascade JCM system in the rotating wave approximation [17,18] is described by the 
Hamiltonian 

H = huia^fa + I,) + (AI, + gh{I^a + I-a^)), (15) 

where and a be the creation and annihilation operators, g be the coupling constant and 
A = h{wQ — w) be the detuning frequency respectively. It is easy to check that both diagonal 
and interaction parts of the Hamiltonian commute with each other. The eigenfunction of this 
Hamiltonian is given by 



li^nit) S [C'l-^\t)\n + l, - > +CS{t)\n,0> +C!l-\t)\n -!,+ >]. 

n=0 



(16) 



We note that the Hamiltonian couples the atom-field states \n — l,+ >, \n,0 > and jn + l, — >, 
where n represents the number of photons of the field. The interaction part of the Hamiltonian 
(15) can also be written in matrix form 



int 



-A gh^/n + 1 
gh\/n + 1 gfi\/n 
gh^/n A 



At resonance (A = 0), the eigenvalues of the Hamiltonian are given by 
and A_ = —gh\j2n + 1 with the corresponding dressed eigenstates 



(17) 

gUsJIn + 1, Ao = 
(18) 



In Eq.(18), the dressed states are constructed by rotating the bare states with the Euler matrix 
T parameterized as 
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where 



"11 = 


cosipcosi^ — cosOsincpsinil) 


"12 = 


cosipsin4> + cos9cos4>sinip 


"13 = 


sinipsinO 


"21 = 


—sinipcos(f) — cosOsirKpcosip 
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—sinipsirKf) + cos9cos(f)cos%l) 


"23 = 


cosipsinO 
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sinOsincf) 
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— sinOcoscj) 
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cos6. 



The evaluation of its various elements is presented in the appendix and here we quote the 
results as follows 
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"31 = ^J 
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The time dependent probability amplitudes of the three levels are given by 
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where Vtn = 5'\/2n + 1. In the following we consider different initial condition of the atom with 

the quantized field in a number state \n >. 

Case-IV: Here we consider the atom is initially polarized in the lower level and the 
combined atom-field state is |n + 1, - >, i.e, C^-^(O) = 0, Co"(0) = 0, C!!+i(0) = 1. Using 
Eqs(20) and (21) the time dependent atomic population of the three levels are given by 



^o(^)l = (2^+1) 



(22a) 
(22b) 

(22c) 



Case-V: At t = when the atom is in the middle level and the combined atom-field state is 
|n, > i.e, C^-i(O) = 0, C^(0) = 1, C!!+i(0) = 0, we find 



C" ^(t) — (2„+i)*^^^^n^, 



C!!+i(t) 



(2n+l) 
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Case- VI: C;-i(0) = 1, Co"(0) = 0, C!!+i(0) = 

If the atom is initially in the upper level and the atom-field state is \n — l,-\- >, i.e, C"~^(0) = 1, 
Q(0) = 0, C!!+^(0) = we obtain the following probabilities 



Cl-\t) = 1 - 4[g^ + ^^cos'nj/2]stn'nj/2, 



C!!+i(t) 



4n(ri+l) - 40 ^/o 



(24a) 
(24b) 

(24c) 



Finally we note that, at resonance, for large value of n the probabilities of case-IV, V and 
VI are identical to case-I, II and III, respectively indicating the validity of the correspondence 
principle. 

IV. Numerical results 

To explore the physical content, we now proceed to analyze the probabilities of the semi- 
classical model and the cascade JCM numerically. 

For the classical field at resonance, the time evolution of the probabilities |C+(t)p (solid 
line), |Co(t)P (dashed line) and |C_(t)p (dotted line) corresponding to Case-I, II and III, 
respectively are shown in Fig 1. We note that for the cases with atom initially in lower and 
upper level, which are displayed in Fig 1(a) and 1(c), respectively the probabilities |C+(t)p 
and |C_(t)p can attain a maximum value equal to unity while |Co(t)P cannot. If we compare 
these two figures the time dependent populations of the lower and upper levels are different by 
a phase lag corresponding to the initial condition of population. This clearly shows that the 
probabilities oscillate between the levels | -I- > and | — > alternatively at a Rabi frequency of 
J^i — On the contrary, the plot of case-II where the atom is initially in the middle level 
depicted in Fig 1(b) shows that the system oscillates with a Rabi frequency of uu — ^ such 
that the probabilities of | + > and | — > states are always equal. When the atom is initially 
in the middle level, the exactly sinusoidal resonant field interacts with the atom in such a way 
that the upper and lower levels are dynamically treated on an equal footing. This dynamically 
symmetrical distribution of population between the upper and lower levels is possible because 
of the classical field. 



For quantized field we consider the time evolution of the probabilities in two different 
situations of initial condition of the field: a) when the field is in a number state and b) the field 
is in a coherent state. 

a) For the cascade JCM, the probabilities of casc-lV, V and VI are plotted in Fig 2 when 
the field is in a number state with n = 1 and g=0.1. In Fig 2(a) we note that for case-IV 
i.e, when the atom is initially in the lower level, the Rabi- frequency of oscillation is = 
However, unhke case-I of the semiclassical model, the probabilities \C'1r^{t)\'^ never becomes 
unity. On the other hand. Fig 2(b) illustrates the probabihties of case-V i.e, when the atom is 
initially in the middle level, where the system oscillates with a Rabi frequency z/^/ — ~ 
once again, in contrast to the corresponding semiclassical situation in case-II, the probabilities 
of the upper and lower level are not equal. The probabilities of case- VI i.e, when the atom 
is initially in the upper level, depicted in Fig 2(c) shows that although it possesses the same 
Rabi frequency uf, but the pattern of oscillation is not out of phase of case-IV. To compare 
with one can look back the semiclassical interaction where we have shown that in case-I the 
pattern of oscillation of upper(lower) level population is precisely identical to the lower(upper) 
level population of case-Ill. 

To understand the implications of such dynamical symmetry breaking qualitatively, var- 
ious bounds on the probabilities are given below : 

Table-I 



Case 


Semiclassical model 


Case 


Cascade JCM 


I 


0< |c_(t)p < 1, 

0<|Co(t)P<l, 
< |C+(t)|2 < 1 


IV 


< |C"(t)p < 1. 

o<|Q(t)p<i, 
< \ci{t)\^ < 1 


II 


< |C_(i)p < 1, 
0<|Co(t)p<l, 
< |C+(t)|2 < 1 


V 


< |C"(i)p < 1, 

o<|Q(t)p<i, 
< |c^(t)p < 1 


III 


< |C_(t)|2 < 1, 
0< |Co(t)P<l, 

< |cv(Oi^ < 1 


VI 


< |C!!(t)|2 < 1, 
< |Co(t)"P < 1, 

< \C"l{t)\'^ < 1 



We note that for the semiclassical model, the symmetric evolution of the probabilities results an 
identical bounds for case-I and III as shown in the Table-I. On quantization of the field mode, 
the bounds corresponding to case-IV and VI are no longer similar although those for case-II 
remains same as case-V. At resonance, for large value of n, Eqs.(22), (23) and (24) of case-IV, 
V and VI are precisely identical to Eqs.(12), (13) and (14) of case-I, II and III, respectively 
and we recover the same bounds of the semiclassical model. 

b) Finally, we consider the atom is interacting with the quantized field mode in a coherent 
state. The coherently averaged probabilities of case-IV, V and VI are given by 



<P+{t) >=Y.Pn Cl-\t) 
n 

<Po{t) >=j:Pn\cmf 



<P-{t) >=EPn Cl+\t) 



(25a) 
(25b) 

(25c) 



where P„ = 



- be the Poisson distribution function and n be the mean photon number. 
For all numerical purpose we choose g — Q.l. We have studied extensively for various values of 



n. The figures are given only for n = 50. Tlie Figs 3-5 display the numerical plots of Eqs.(25) 
where the collapse and revival of the Rabi oscillation is clearly evident. For low n and when the 
atom is in the middle level the symmetrical values of population of upper and lower levels are 
not observed until n is very high as given in the figures. However, even if n = 50, the numerical 
values of the time dependent populations of the upper and lower levels are not exactly equal 
although very close and becomes exactly equal in the limit n — > oo. We further note that, 
if the atom is initially polarized either in the upper or in the lower level, it exhibits similar 
population oscillation, which is different from the case if it is initially polarized in the middle 
level. The reproduction of this result analogous to the semiclassical model shows the proximity 
of the coherent state with large n to the classical field. 

When the field is quantized, the population oscillation depends on the occupation number, 
n of the field state, for example, cos {g\/1n + It). For a statistical distribution of field state, 
the spontaneous factor 1 plays a dominant role when n is low. For an initial number state of 
the field when n is slightly higher than 1, the upper and lower levels of the atom are not treated 
dynamically on an equal footing even when the atom is initially in the middle level. This fine 
graining of the quantized distribution of photons over the number states {jri >} generates a 
complex interference between individual Rabi oscillations corresponding to each n and plays 
a role until when n is very large compared to 1 and effectively acts as a classical field and 
thereby the semiclassical situation is satisfied. Note that for an initial vacuum field, i.e, n = 
for the number state and n = for the coherent state, with the atom initially in the middle 
level, it can not go to the upper level at all and the population will oscillates between the lower 
and middle levels with Rabi frequency Oq- This asymmetry is still present when the field is 
in a coherent state with a Poissonian photon distribution with low average photon number, n, 
which is generally not symmetric around an n. A Poisson distribution is almost symmetric, a 
Gaussian, around a n if n is very large which is the case of a classical field. In that situation the 
upper and lower levels of the atom are treated dynamically on an equal footing and maintains 
the symmetrical distribution of population in upper and lower levels. 

V. Conclusion 

We conclude by recapitulating the essential content of our investigation. At the outset 
we have sculpted the semiclassical model by choosing the spin-one representation of SU{2) 
group and have calculated the transition probabilities of the three levels. It is shown that at 
resonance, if the atom is initially polarized in the lower or in the upper level, the various atomic 
populations oscillate quite differently when it is initially populated in the middle level. When 
the atom is initially populated in the middle level, the classical field interacts in such a way that 
the populations of the upper and lower levels are always equal. This dynamically symmetrical 
populations of the upper and lower levels are destroyed due to the quantization of the field. To 
show this quantum behaviour a dressed-atom approach is presented to solve the cascade JCM. 
Finally we discuss the restoration of the symmetry taking the quantized field in a coherent state 
with large average photon number, a closest state to the classical state. Although the collapse 
and revival and some other nonclassical features are well studied in the context of two-level 
systems, the above dynamical breaking of symmetry due to the quantization of the field has 
no two-level analog. We hope that this dynamical behaviour in the cascade three level system 
should show its signature on the time dependent profile of the second order coherence of the 
quantized field which will be discussed elsewhere. The dressed-atom approach developed here 
may also find its application in the V and A type three level systems where the nature of the 
symmetry should be different from the cascade system. 



APPENDIX 



At resonance, the interaction part of the Hamiltonian of the three level system is given by 



int 



9\/nTl 
g\/n + 1 g^Jn 
gy/n 

= g\J1n + 1, Ao = and A_ 



(A.l) 



where the eigenvalues are A+ = g\/1n^\^ Aq = and A_ = —g\/2n + 1. The Euler matrix 
T, diagonahzes the Hamiltonian as Ho — THintT~^, is given by Eq.(19). Using the trick 
{Hint — XjI){Xj} — 0, where {Xj} be the column matrix of T"^, corresponding to the eigenvalue 
X+ we have 



-gV^n + 1 gVn + 1 
gy/n + 1 -gy/2n + 1 g ^/n 
S'V^ — g'V'SjT' + 1 





an 




"12 







= 0. 



These linear equations readily yield 

\/2n+l 



q;i2 



"13, "12 



V2re+1 
Vn+T 



an and an 



/n+l 



ai3. 



(A.2) 



(A.3) 



Using the normalization condition 

«11 + «12 + "13 

we get an = a/^J , "12 = 



2n+l 



^ and ai3 



4n+2 V2 V 4n+2 

eigenvalues Aq and A_ we can obtain other elements of T, namely, 



(A.4) 

. Similarly, corresponding to the 



•^21 = --/iSt ' «22 = and a23 = \j 



n+l 
2n+l' 



«3i = "32 = -71 and a33 = 

One can now easily read off the Euler's angles 



4n+2' 



= a/S' ^''^'^ ^ \/S = a/52- 



(A.5) 
(A.6) 

(A.7) 
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[Fig.l] : The time evolution of the probabiUties of the semiclassical model corresponding 
to case-I, II and III. The symmetric pattern of evolution is evident from Fig. la and 
Fig.lc which are in opposite phase. 
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[Fig.2] : The time evolution of the probabilities of the cascade JCM corresponding to 
case-IV, V and VI. In Fig.2a and Fig.2c depict that the symmetry exhibited by the 
semiclassical model is spoiled on quantization of the field mode. 
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[Fig.3] : The collapse and revival are shown for case-IV when the atom is initially in the 
lower level. The time dependent profiles of the upper and lower level populations are 
similar in Fig 3(a) and 3(c). 



Fig. 4a Fig. 4b Fig. 4c 
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[Fig 4] : The collapse and revival for case-V are displayed when the atom is initially in 
the middle level. The time dependent patterns of the upper and lower level populations 
are similar in Fig 4(a) and 4(c) as in the semiclassical cases of Fig 1(b). 
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[Fig5] : The collapse and revival for case- VI are shown when the atom is initially polarized 
in the upper level. Similar to Fig 1(a) and Fig 1(c) of the semiclassical model, Fig 3 and 
Fig 5 are closely alike. 



